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We report the linear optical absorption spectra of aluminum clusters Aln (n=2-5), using system- 
atic large-scale all-electron correlated calculations. Several low- lying isomers of each cluster were 
considered, and their geometries were optimized at the coupled-cluster singles doubles (CCSD) level 
of theory. With these optimized ground-state geometries, excited states of different clusters were 
computed using the multi-reference singles-doubles configuration-interaction (MRSDCI) approach, 
which includes electron correlation effects at a sophisticated level. These CI wave functions were 
used to compute the transition dipole matrix elements connecting the ground and various excited 
states of different clusters, eventually leading to their linear absorption spectra. The convergence of 
our results with respect to the basis sets, and the size of the CI expansion, was carefully examined. 
The contribution of configurations to many body wavefunction of various excited states suggests 
that the excitations involved are collective, and plasmonic in nature. 
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I. INTRODUCTION 

Metal clusters are promising candidates in the era of 
nanotechnology. The reason behind growing interest in 
this field lies in interesting properties, applicability of 
very simple theoretical models to describe their proper- 
ties, and a vast variety of potential technological appli- 
cations [1-5]. 

Various jellium models have successfully described 
electronic structures of alkali metal clusters, because al- 
kali metals have free valence electrons. This beautifully 
explains the magic abundance of certain clusters. How- 
ever, in case of aluminum clusters, the experimental re- 
sults often provide conflicting evidence about the size at 
which the jellium model would work. The theoretical 
explanation also depends on the valency of aluminum 
atoms considered. Since s - p orbital energy separation 
in aluminum atom is 4.99 eV, and it decreases with the 
cluster size, the valency should be changed from one to 
three [6]. Perturbed jellium model, which takes orbital 
anisotropy into account, has successfully explained the 
mass abundance of aluminum clusters [7, 8]. 

Shell structure and s - p hybridization in anionic alu- 
minum clusters were probed using photoelectron spec- 
troscopy by Grantefor and Eberhardt [9], and Li et al. 
[10]. Evolution of electronic structure and other prop- 
erties of aluminum clusters has been studied in many 
reports [6, 7, 10-20]. Structural properties of aluminum 
clusters were studied using density functional theory by 
Rao and Jena [6]. An all electron and model core poten- 
tial study were carried out by Martinez et al. [19]. Upton 
performed chemisorption calculations on aluminum clus- 
ters and reported that Ale will be the smallest cluster 
that will absorb H2 . DFT alongwith molecular dynamics 
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were used to study electronic and structural properties 
of aluminum clusters [16]. 

The optical absorption in alkali metal clusters has been 
studied extensively at various levels of theory [21-24]. 
On the contrary, although there is enormous informa- 
tion available about ground state structures, very lit- 
tle information is available about excited states of alu- 
minum clusters. Optical absorption in several aluminum 
clusters corresponding to the minimum energy configu- 
rations has been studied by Deshpande et al. using time- 
dependent density functional theory (TDDFT) [25]. Xie 
et al. presented TDDFT optical absorption spectra of 
various caged icosahedral aluminum clusters [26]. How- 
ever, no theoretical or experimental study has been done 
on optical absorption in various low lying isomers of alu- 
minum clusters. The distinction of different isomers of a 
cluster has to be made using an experimental or theoret- 
ical technique in which the properties are size, as well as 
shape, dependent. Conventional mass spectrometry only 
distinguishes clusters according to the masses. Hence, 
our technique can be coupled with the experimental mea- 
surements of optical absorption to distinguish between 
different isomers of a cluster. We have recently reported 
results of such calculation on small boron clusters [27]. 

In this paper, we present results of systematic calcula- 
tions of optical absorption in various low lying isomers of 
small aluminum clusters using ah initio large-scale con- 
figuration interaction method. The nature of optical ex- 
citations involved in absorption have also been presented 
by analyzing the wave functions of the excited states. 

Remainder of the paper is organized as follows. Next 
section discusses theoretical and computational details 
of the calculations, followed by section HI, in which re- 
sults are presented and discussed. Conclusions and fu- 
ture directions are presented in section IV. Paper ends 
with an appendix presenting detailed information about 
wavefunctions of excited states contributing to various 
photoabsorption peaks. 
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II. THEORETICAL AND COMPUTATIONAL 
DETAILS 

The geometry of various isomers were optimized us- 
ing the size-consistent coupled-cluster singles-doubles 
(CCSD) method, as implemented in the Gaussian 09 
package [28]. A basis set of 6-311++G(2d,2p) was used 
which was included in the Gaussian 09 package itself. 
This basis set is optimized for the ground state calcula- 
tions. Since an even numbered electron system can have 
singlet, triplet, or higher spin multiplicity, we repeated 
the optimization for singlet and triplet systems to look for 
true ground state geometry. Similarly, for odd numbered 
electron systems, doublet and quartet multiplicities were 
considered in the geometry optimization. To initiate the 
optimization, raw geometries, reported by Rao and Jena, 
based on density functional method were used [6]. Fig- 
ure 1 shows the final optimized geometries of the isomers 
studied in this paper. 

Using these optimized geometries, correlated calcula- 
tions were performed using multireference singles doubles 
configuration interaction (MRSDCI) method for both 
ground state and excited states [29]. This method con- 
siders a large number of singly- and doubly- substituted 
configurations from a large number of reference configu- 
rations, and, is well suited for both ground and excited 
states calculations. It takes into account the electron cor- 
relation which is inadequately represented in single refer- 
ence ab initio methods. These ground- and excited-state 
wavefunctions are further used to calculate the transi- 
tion dipole moment matrix elements, which in turn, are 
utilized to compute linear optical absorption spectrum 
assuming a Lorentzian line shape. 

Various wave functions of the excited states contribut- 
ing to the peaks in the spectrum were analyzed, and even 
bigger MRSDCI calculations were performed by includ- 
ing more references, if needed. The criteria of choosing 
a reference in the calculation was based upon the mag- 
nitude of the corresponding coefficients in the CI wave 
function of the excites states contributing to a peak in the 
spectrum. This process was repeated until the spectrum 
converges within acceptable tolerance and all the config- 
urations which contribute to various excites states were 
included. The typical number of reference configurations 
and total number of configurations considered in the cal- 
culations of various isomers are given in Table I. We have 
extensively used such approach in performing large-scale 
correlated calculations of linear optical absorption spec- 
tra of conjugated polymers [30-33]. and atomic clusters 
[27]. 

Since configuration interaction method is computa- 
tionally very expensive, mainly, because the number of 
determinants to be considered increases exponentially 
with the number of electrons, and the number of molec- 
ular orbitals. Calculations on bigger clusters are pro- 
hibitive under such circumstances, and are very time 
consuming even for the clusters considered here. There- 
fore, whenever possible, efforts were made to minimize 



the computational time by using systems properties and 
some approximations. Point group symmetries (D2/1, and 
its subgroups) were taken into account, thereby making 
calculations for each symmetry subspace independent of 
each other. The core of the aluminum atom were frozen 
from excitations, keeping only three valence electrons ac- 
tive. Also an upper limit on the number of virtual or- 
bitals was imposed, to restrict very high energy excita- 
tions. The effect of these approximations on the com- 
puted photoabsorption spectra has been studied carefully 
and is presented in the next section. 



III. RESULTS AND DISCUSSION 

In this section, we have presented a systematic study 
of the convergence of our results and various approxima- 
tions used. In the latter part, we discuss the results of 
our calculations of various clusters. 



A. Convergence of calculations 

In this section we discuss the convergence of photoab- 
sorption calculations with respect to the choice of the 
basis set, and the size of the active orbital space. 



i. Choice of basis set 

It is well known that, electronic structure calculations 
are basis set dependent. Several basis sets are optimized 
for specific purposes, such as ground state calculations, 
excited states and covalent systems etc. We have re- 
ported a systematic basis set dependence of photoab- 
sorption of boron cluster [27]. Similarly, here we have 
checked the dependence of photoabsorption spectrum of 
aluminum dimer on basis sets used [34, 35] as shown in 
Fig. 2. The 6-311 type Gaussian contracted basis sets are 
known to be good for ground state calculations. The cor- 
relation consistent (CC) basis sets, namely, CC-polarized 
valence double-zeta and CC-polarized valence triple zeta 
(cc-pVTZ) give a good description of excited states of 
systems. The latter is found to be more sophisticated in 
describing the high energy excitations, which were also 
confirmed using results of an independent TDDFT cal- 
culation [36]. So we have used cc-pVTZ basis set for the 
optical absorption calculations. 



2. Orbital truncation scheme 

With respect to the total number of orbitals N in the 
system, the computational time in configuration interac- 
tion calculations scales as ^ N^. Therefore such calcu- 
lations become intractable for moderately sized systems, 
such as those considered here. So, in order to ease those 
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Figure 1: (Color online) Geometry optimized structures of aluminum clusters with point group symmetry and the 
electronic ground state at the CCSD level. All numbers are in A unit. 
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Figure 2: (Color online) Optical absorption in AI2 calcu- 
lated using various Gaussian contracted basis sets. 



calculations, the lowest lying molecular orbitals are con- 
strained to be doubly occupied in all the configurations, 
implying that no optical excitation can occur from those 
orbitals. It reduces the size of the CI Hamiltonian ma- 
trix drastically. In fact, this approach is recommended 
in quantum chemical calculations, because the basis sets 
used are not optimized to incorporate the correlations in 
core electrons [37]. The effect of this approximation on 
the spectrum is as shown in Fig. 3. Since, calculations 
with all electrons in active orbitals were unfeasible, we 
have frozen occupied orbitals upto -4 Hartree of energy 
for the purpose of demonstration. The effect of freezing 
the core is negligibly small in the low energy regime, but 
shows disagreement in the higher energy range. How- 
ever, for very high energy excitations, photodissociation 
may occur, hence absorption spectra at those energies 



Figure 3: (Color online) The effect of freezing the core 
orbitals of aluminum atoms on optical absorption spec- 
trum of AI2. It renders little effect on optical absorp- 
tion spectrum, with significant reduction in the compu- 
tational cost. 



will cease to have meaning. Thus, the advantage of freez- 
ing the core subdues this issue. Therefore, in all the cal- 
culations, we have frozen the chemical core from optical 
excitations. 

Not only occupied, but high energy virtual (unoccu- 
pied) orbitals can also be removed from the calculations 
to make them tractable. In this case the high lying or- 
bitals are constrained to be unoccupied in all the config- 
urations. This move is justifiable, because it is unlikely 
that electrons would prefer partial filling of high energy 
orbitals in an attempt to avoid other electrons. However, 
this will only be applicable if the orbitals are sufficiently 
high in energy. Fig. 4 shows the effect of removing or- 
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Figure 4: (Color online) The effect of the number of ac- 
tive orbitals (N^ct) on the optical absorption spectrum 
of AI2. Until Nac^=46, the optical spectrum does not 
exhibit any significant change. It corresponds to 1.0 
Hartree {^27.2 eV) virtual orbital energy. 



bitals having more than the specified energy. From the 
figure it is clear that photoabsorption spectra exhibits no 
difference at all upto 1 Hartree cutoff on virtual orbitals. 
Below 0.8 Ha cutoff, the spectra start deviating from each 
other. Hence, we have ignored the virtual orbitals having 
energy more than 1 Ha. 



3. Size of the CI expansion 

In the multireference configuration interaction method, 
the size of the CI Hamiltonian increases exponentially 
with the number of molecular orbitals in the system. 
Also, accurate correlated results can only be obtained 
if sufficient number of reference configurations are in- 
cluded in the calculations. In our calculations, we have 
included those configurations which are dominant in the 
wave functions of excited states for a given absorption 
peak. Also, for ground state calculations, we included 
configurations until the total energy converges within a 
pre-defined tolerance. Table I shows the average num- 
ber of reference configurations and average number of 
total configurations involved in the CI calculations of 
various isomers. For a given isomer, the average is calcu- 
lated across different irreducible representations needed 
in these symmetry adapted calculations of the ground 
and various excited states. For the simplest cluster, the 
total configurations are about half a million and for the 
biggest cluster considered here, it is around four million 
for each symmetry subspace of AI5. The superiority of 
our calculations can also be judged from the correlation 
energy defined here (c/. Table I), which is the difference 
in the total energy of a system at the MRSDCI level and 
the Hartree- Fock level. The correlation energy per atom 
seems to be quite high for all the clusters, making our 
calculations stand out among other electronic structure 



Figure 5: (Color online) Molecular orbitals of aluminum 
dimer. H and L stands for HOMO and LUMO respec- 
tively, and Hi and H2 are singly occupied degenerate 
molecular orbitals. 



calculations, especially single reference DFT based cal- 
culations. 



B. Calculated photoabsorption spectra of various 
clusters 

In this section, we describe the photoabsorption spec- 
tra of various isomers of the aluminum clusters studied. 
Graphical presentation of molecular orbitals involved are 
also given in each subsection below. 

L Ah 

Aluminum dimer is the most widely studied cluster of 
aluminum. It has now been confirmed that the AI2 (c/. 
Fig. 1(a)) has ^H^ ground state. The bond length ob- 
tained using geometry optimization at CCSD level was 
2.72 A with Doo/i point group symmetry. This is in very 
good agreement with available data, such as Martinez et 
al. obtained 2.73 A as dimer length using all electron cal- 
culations [19], 2.71 A [17] and 2.75 A [38] as bond lengths 
using DFT and configuration interaction methods, and 
2.86 A obtained using DFT with generalized gradient 
approximation [6]. The experimental bond length of alu- 
minum dimer is 2.70 A [39]. Another metastable state of 
the dimer exists with electronic state, and 2.48 A in 
bond length. 

The many-particle wave function of AI2 consists of two 
degenerate singly occupied molecular orbitals (to be de- 
noted by Hi and i^2, henceforth), because it is a spin 
triplet system. Similarly, the configurations involving ex- 
citations from occupied molecular orbitals to the unoc- 
cupied orbitals, form excited state wave functions. The 
computed photoabsorption spectra of AI2, as shown in 
Fig. 6, is characterized by weaker absorption at lower en- 
ergies and couple of intense peaks at higher energies. The 
many-particle wave functions of excited states contribut- 
ing to the peaks are presented in Table A.l. The spec- 
trum starts with a small absorption peak at around 2 eV, 
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Table I: The average number of reference configurations (N^g/), and average number of total configurations total) 
involved in MRSDCI calculations, ground state (GS) energies (in Hartree) at the MRSDCI level, relative energies and 
correlation energies (in eV) of various isomers of aluminum clusters. 



Cluster 


Isomer 




^total 


GS energy 
(Ha) 


Relative 
energy (eV) 


Correlation energy^ 
per atom(eV) 




T ,1 n T* 


40 


445716 


'-too .tyyjKJO^o ± 


0.0 


1.69 


AI3 


Equilateral triangular 


40 


1917948 


-725.9053663 


0.00 


2.38 




Isosceles triangular 


22 


1786700 


-725.8748996 


0.83 


2.36 




Linear 


18 


1627016 


-725.8370397 


1.85 


2.16 


AI4 


Rhombus 


13 


3460368 


-967.8665897 


0.00 


1.82 




Square 


21 


1940116 


-967.8258673 


1.11 


1.80 


AI5 


Pentagonal 


7 


3569914 


-1209.8114803 


0.00 


1.73 




Pyramidal 


8 


3825182 


-1209.7836568 


0.76 


1.77 



^ The difference in Hartree-Fock energy and MRSDCI correlated energy of the ground state. 



characterized by ^ ^+1 and light polarized along the 
direction of axis of the dimer. It is followed by a couple of 
small intensity peaks, until a dominant absorption is seen 
at 5 eV. This is characterized by i^i ^ L + 3. Another 
dominant peak is observed at 8 eV having H — 2 ^ L 
as dominant configuration, with absorption due to light 
polarized perpendicular to the axis of the dimer. 

If we compare this spectrum with the one ob- 
tained from time-dependent local density approximation 
(TDLDA) method [25], we get a fair agreement - al- 
though there is disagreement with respect to the intensity 
and number of peaks. Unlike our calculations, the num- 
ber of peaks is much more in TDLDA results and the 
spectrum is almost continuous. Two major peaks at 5 eV 
and 8 eV are in complete agreement with the previously 
reported results of Deshpande et al. [25] . Minor peaks at 
3.2 eV and 6.3 eV are also observed in the TDLDA spec- 
trum of dimer. The latter one is, contrary to TDLDA 
spectrum, found to be a small peak in our calculations. 



2. Ak 

Among the possible isomers of aluminum cluster AI3, 
the equilateral triangular isomer is found to be the most 
stable. We have considered three isomers of AI3, namely, 
equilateral triangle, isosceles triangle, and a linear chain. 
The most stable isomer has D^h point group symmetry, 
and '^Ai electronic state. The optimized bond length 2.57 
A, is in good agreement with reported theoretical values 
2.61 A [6], 2.62 A [7], 2.56 A [19], and 2.52 A [12, 13]. The 
doublet ground state is also confirmed with the results of 
magnetic deflection experiments [15]. 

The next isomer, which lies 0.83 eV high in energy, 
is isosceles triangular isomer. The optimized geometry 
has 2.59 A, 2.59 A and 2.99 A as sides of triangle, with a 
quartet ground state (^^2). Our results are in agreement 
with other theoretical results [7, 16, 19]. 




Figure 6: (Color online) The linear optical absorption 
spectrum of AI2, calculated using the MRSDGI approach. 
The peaks corresponding to the light polarized along 
the molecular axis are labeled with the subscript ||, 
while those polarized perpendicular to it are denoted by 
the subscript ±. For plotting the spectrum, a uniform 
linewidth of 0.1 eV was used. 



Linear AI3 isomer again with quartet multiplicity is the 
next low- lying isomer. The optimized bond length is 2.62 
A. This is in good agreement with few available reports 
[12, 16, 19]. 

The photoabsorption spectra of these isomers are pre- 
sented in Figs. 7, 9 and 11. The corresponding many 
body wave functions of excited states corresponding to 
various peaks are presented in Table A. 2, A. 3 and A. 4 
respectively. In the equilateral triangular isomer, most 
of the intensity is concentrated at higher energies. The 
same is true for the isosceles triangular isomer. How- 
ever, the spectrum of isosceles triangular isomer appears 
slightly red shifted with respect to the equilateral coun- 
terpart. Along with this shift, there appears a split pair 
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Figure 7: (Color online) The linear optical absorption 
spectrum of AI3 equilateral triangle isomer, calculated 
using the MRSDCI approach. The peaks corresponding 
to the light polarized along the molecular plane are la- 
beled with the subscript ||, while those polarized perpen- 
dicular to it are denoted by the subscript ±. For plotting 
the spectrum, a uniform linewidth of 0.1 eV was used. 



Figure 9: (Color online) The linear optical absorption 
spectrum of AI3 isosceles triangle isomer, calculated us- 
ing the MRSDCI approach. All peaks labeled above cor- 
respond to the light polarized along the molecular plane. 
For plotting the spectrum, a uniform linewidth of 0.1 eV 
was used. 



L+3 L+2 L+1 L H 
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Figure 8: (Color online) Molecular orbit als of equilateral 
triangular aluminum trimer. H and L stands for HOMO 
and LUMO respectively. 



of peaks at 5.8 eV. This splitting of oscillator strengths 
might be due to distortion accompanied by symmetry 
breaking. The absorption spectrum of linear isomer is 
altogether different with bulk of the oscillator strength 
carried by peaks in the range 4-5 eV, and, due to the 
polarization of light absorbed parallel to the axis of the 
trimer. 

The optical absorption spectrum of equilateral trian- 
gular isomer consists of very feeble low energy peaks at 
3.5 eV, 5.6 eV and 5.8 eV characterized by — 3 ^ L + 5, 
a double excitation i^ — 2^L + 5;i^ — l^L + 5, and 
H — 3 L + 2 respectively. The latter peak is due to the 
light polarized perpendicular to the plane of the isomer. 
It is followed by an intense peak at around 6.5 eV with 
dominant contribution from H ^ L -\- 6 and H ^ L -\- 4 
configurations. A semi-major peak is observed at 7.5 eV 
characterized mainly due to double excitations. 

Two major peaks at 6.5 eV and 7.5 eV, obtained in our 
calculations are also found in the spectrum of TDLDA 



calculations, with the difference that the latter does not 
have a smaller intensity in TDLDA [25]. Other major 
peaks obtained by Deshpande et at. [25] in the spectrum 
of aluminum trimer are not observed, or have very small 
intensity in our results. 

As compared to the equilateral triangle spectra the 
isosceles triangular isomer exhibit several small intensity 
peaks (c/. Fig. 9) in the low energy regime. The ma- 
jority of contribution to peaks of this spectrum comes 
from longitudinally polarized transitions, with negligible 
contribution from transverse polarized light. The spec- 
trum starts with a feeble peak at 2.4 eV with contribution 
from doubly-excited configuration H^L-\-l;H — 2^ 
L + 2. One of the dominant contribution to the oscil- 
lator strength comes from two closely- lying peaks at 5.8 
eV. The wave functions of excited states corresponding 
to this peak show a strong mixing of doubly-excited con- 
figurations, such asJL — 3 L + l;i7 — 2 L and 
H-2^L^1;H-4^L. The peak at 6.7 eV shows 
absorption mainly due to JL ^ L + 10. 

Linear trimer of aluminum cluster also shows low ac- 
tivity in the low energy range. Very feeble peaks are 
observed at 1.2 eV and 2.3 eV, both characterized by 
H — 3 H — 2. This configuration also contributes to 
the semi-major peak at 4 eV along with H — 4 ^ H. Two 
closely lying peaks at 4.3 eV and 4.6 eV carry the bulk of 
the oscillator strength. Major contribution to the former 
comes from H — 1 ^ L + 2 along with H — 3 ^ H — 2 
being dominant in both the peaks. Again, as expected, 
the absorption due to light polarized along the trimer 
contributes substantially to the spectrum. 
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Figure 10: (Color online) Molecular orbitals of isosceles 
triangular aluminum trimer. H and L stands for HOMO 
and LUMO respectively, and i^i, i^2, and H^, are singly 
occupied molecular orbitals. 




Figure 11: (Color online) The linear optical absorption 
spectrum of AI3 linear isomer, calculated using the MRS- 
DCI approach. The peaks corresponding to the light 
polarized along the molecular axis are labeled with the 
subscript ||, while those polarized perpendicular to it are 
denoted by the subscript ±. For plotting the spectrum, 
a uniform linewidth of 0.1 eV was used. 




H-2 



Figure 12: (Color online) Molecular orbitals of linear alu- 
minum trimer. H and L stands for HOMO and LUMO 
respectively, and i^i, i^2, and H^, are singly occupied 
molecular orbitals. 




Figure 13: (Color online) The linear optical absorption 
spectrum of rhombus AI4, calculated using the MRSDCI 
approach. The peaks corresponding to the light polarized 
along the molecular axis are labeled with the subscript 
II, while those polarized perpendicular to it are denoted 
by the subscript ±. For plotting the spectrum, a uniform 
linewidth of 0.1 eV was used. 



3. Ah 

Tetramer of aluminum cluster has many low lying iso- 
mers due to its flat potential energy curves. Among them, 
rhombus structure is the most stable with ^B2g electronic 
ground state. Our optimized bond length for rhombus 
structure is 2.50 A and 63.8° as the acute angle. This 
is to be compared with corresponding reported values of 
2.56 A and 69.3° reported by Martinez et al [19], 2.51 A 
and 56.5° computed by Jones [17], 2.55 A and 67.6° ob- 
tained by Schultz et al [13]. We note that bond lengths 
are in good agreement but bond angles appear to vary a 
bit. 

The other isomer studied here is an almost square 
shaped tetramer with optimized bond length of 2.69 A. 
The electronic ground state of this D/^h symmetric clus- 
ter is ^^sn- This optimized geometry is in accord with 
2.69 A reported by Martinez et al. [19], however it is 
somewhat bigger than 2.57 A calculated by Yang et al. 
[12] and 2.61 A obtained by Jones [16]. 

For planar clusters, like rhombus and square shaped 
AI4, two types of optical absorptions are possible: (a) 
Longitudinal - those polarized in the plane of the cluster, 
and (b) transverse - the ones polarized perpendicular to 
that plane. The onset of optical absorption in rhombus 
isomer occurs at around 1 eV with transversely polarized 
absorption characterized by i^i ^ L + 1. It is followed 
by an in-plane polarized absorption peak at 2.3 eV with 
dominant contribution from H — 2 Hi. Several closely 
lying peaks are observed in a small energy range of 4.5 - 
8 eV. Peaks split from each other are seen in this range 
confirming that after shell closure, in perturbed droplet 
model, Jahn Teller distortion causes symmetry breaking 
usually associated with split absorption peaks. The most 
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Figure 14: (Color online) Molecular orbitals of rhombus- 
shaped aluminum tetramer. H and L stands for HOMO 
and LUMO respectively, and Hi and H2 are singly occu- 
pied molecular orbitals. 
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Figure 15: (Color online) Molecular orbitals of square- 
shaped aluminum tetramer. H and L stands for HOMO 
and LUMO respectively, and Hi and H2 are singly occu- 
pied molecular orbitals. 




Figure 16: (Color online) The linear optical absorption 
spectrum of square AI4, calculated using the MRSDCI 
approach. The peaks corresponding to the light polarized 
along the molecular plane are labeled with the subscript 
II, while those polarized perpendicular to it are denoted 
by the subscript ±. For plotting the spectrum, a uniform 
linewidth of 0.1 eV was used. 



The TDLDA spectrum [25] of aluminum rhombus 
tetramer is in very good agreement with our results. 
With a few exceptions of absence of small peaks in very 
low energy range, the spectrum matches very well with 
our spectrum. Peaks labeled HI to XII in our calculated 
spectrum are also observed with TDLDA calculations, 
however, the relative intensity is not always found to be 
consistent. For example, the strongest absorption peak 
of TDLDA calculations, is at around 7.9 eV while we have 
observed a peak at this energy but with small intensity. 
Also the highest absorption peak in our calculations is 
at 5.5 eV, and TDLDA also reports a strong peak at the 
same energy [25]. 



intense peak is observed at 5.5 eV characterized hy H — 
3 ^ L + 4. 

The absorption spectrum of square shaped isomer be- 
gins with a couple of low longitudinally polarized ab- 
sorption peaks at 2.1 eV and 2.7 eV characterized by 
H — 1 ^ L and H2 ^ L -\- 1 respectively. The peak 
at 4.2 and 4.9 eV have H -2 ^ L and Hi ^ L + 2 as re- 
spective dominant configurations. A major peak at 5.85 
eV is observed with absorption due to longitudinal po- 
larization having H — 2 ^ L -\-2 and a double excitation 
Hi L -\- 2] H — 2 L + 2 as dominant configura- 
tions. These configurations also make dominant contri- 
bution to the peak at 6.5 eV. This peak along with one 
at 6.9 eV are two equally and most intense peaks of the 
spectrum. The latter has additional contribution from 
Hi L -\- 1; H — 2 ^ L. A shoulder peak is observed at 
7.2 eV. 



The lowest lying pentagon isomer of aluminum has C2V 
symmetry and has an electronic ground state of^Ai. The 
bond lengths are as shown in Fig. 1(g). These are slightly 
bigger than those obtained by Rao and Jena [6] and Yang 
et al. [12] using the DFT approach. Many other reports 
have confirmed that the planar pentagon is the most sta- 
ble isomer of Al^. 

The other optimized structure of pentamer is perfect 
pyramid with C/iy symmetry and '^Ai electronic ground 
state. This lies 0.76 eV above the global minimum struc- 
ture. This is the only three dimensional structure studied 
in this paper for optical absorption. The optimized ge- 
ometry is consistent with those reported earlier by Jones 
[16]. However, it should be noted that there exists many 
more similar or slightly distorted structure lying equally 
close the the global minimum. 

The optical absorption spectrum of pentagonal AI5 has 




Figure 17: (Color online) The linear optical absorption 
spectrum of pentagonal AI5, calculated using the MRS- 
DCI approach. The peaks corresponding to the light 
polarized along the molecular axis are labeled with the 
subscript ||, while those polarized perpendicular to it are 
denoted by the subscript ±. For plotting the spectrum, 
a uniform linewidth of 0.1 eV was used. 
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Figure 18: (Color online) Molecular orbitals of pentag- 
onal aluminum pentamer. H and L stands for HOMO 
and LUMO respectively. 



few low energy peaks followed by major absorption at 
4.4 eV. It has dominant contribution from H — 1 ^ L -\- 
5 configuration. Pentagonal isomer shows more optical 
absorption in the high energy range, with peaks within 
regular intervals of energy. 

Few feeble peaks occur in the low energy range in the 
optical absorption of pyramidal isomer. The major ab- 
sorption peak at 4.2 eV is slightly red-shifted as com- 
pared to the pentagonal counterpart. It is characterized 
by H — 3 ^ L -\- 2. A peak at 6 eV is seen in this 
absorption spectrum having dominant contribution from 
H 1/ + 13, which is missing in the spectrum of pen- 
tagon. These differences can lead to identification of iso- 
mers produced experimentally. 
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Figure 19: (Color online) Molecular orbitals of pyramidal 
aluminum pentamer. H and L stands for HOMO and 
LUMO respectively. 
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Figure 20: (Color online) The linear optical absorption 
spectrum of pyramidal AI5, calculated using the MRS- 
DCI approach. The peaks corresponding to the light po- 
larized along the Cartesian axes are labeled accordingly. 
For plotting the spectrum, a uniform linewidth of 0.1 eV 
was used. 



The TDLDA calculated spectrum [25] of pentagonal 
isomer is also in very good agreement with our results. 
A small peak at 2.4 eV is observed in both the spectra, 
followed by peaks at 3.9 eV, 4.2 eV and 4.4 eV. These 
three peaks are also observed in TDLDA results with 
a little bit of broadening. Again, the peak at 5.4 eV 
matches with each other calculated from both the ap- 
proaches. Peak found at 6.7 eV is also observed using 
TDLDA calculation. In the range of spectrum studied in 
our calculations, the strongest peak position and inten- 
sity is in accord with that of TDLDA counterpart [25] . 



IV. CONCLUSIONS AND OUTLOOK 

In this study, we have presented large-scale all-electron 
correlated calculations of optical absorption spectra of 
several low- lying isomers of aluminum clusters AU, 
(n=2-5). Both ground and excited state calculations 
were performed at MRSDCI level, which take electron 
correlations into account at a sophisticated level. We 
have analyzed the nature of low-lying excited states. Iso- 
mers of a given cluster show a distinct signature spec- 
trum, indicating a strong-structure property relationship. 
This fact can be used in experiments to distinguish be- 
tween different isomers of a cluster. Owing to sophis- 
tication of our calculations, our results can be used for 
benchmarking of the absorption spectra. The optical ex- 
citations involved are found to be collective type, and 
plasmonic in nature. This can also be used as another 
criterion for evolution of metallic character in metal atom 
clusters. 
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Appendix A: Excited State CI Wavef unctions, 
Energies and Oscillator Strengths 

In the following tables, we present the excitation ener- 
gies (with respect to the ground state), the many body 
wavefunction and the oscillator strengths of the excited 
states corresponding to the peaks in the photoabsorption 
spectra of various isomers listed in Fig. 1, and discussed 
in section III. 
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Table A.l: Excitation energies {E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI2 {cf. Fig. 6), along with the oscillator strength (/12) of the transitions. 
Longitudinal and transverse polarization corresponds to the absorption due to light polarized along and 
perpendicular to the molecular axis respectively. In the wave function, the bracketed numbers are the CI coefficients 
of a given electronic configuration. Symbols i^i,i^2 denote SOMOs discussed earlier, and and L, denote HOMO 
and LUMO orbitals respectively. HF denotes the Hartree-Fock configuration. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 










\Hl, Hi) (0.9096) 

\H -l^ Hi-H2 L)(0.1139) 

\H -2^ E,H -l^ L + 2)(0.0889) 


I 


1.96 


0.1027 


longitudinal 


\H2 ^ L + 1)(0.8120) 
\H -l^ ifi)(0.3685) 


II 


3.17 


0.1249 


longitudinal 


|if-l^i7i)(0.6172) 
\Hi ^ L + 3) (0.4068) 
\Hi ^ L; if- l^L) (0.3190) 


III 


4.47 


0.5149 


transverse 


\H2 ^ L + 4) (0.8313) 
\H2 ^ L + 6) (0.2024) 


IV 


4.99 


5.4531 


longitudinal 


\Hi ^ L + 3) (0.7353) 
if- l^ifi) (0.4104) 


V 


6.31 


0.2554 


transverse 


\H2 ^ L + 6) (0.4683) 

if- l^L + 1) (0.3894) 

if - 1 ^ L; ifs ^ iv + 2) (0.3886) 


VI 


7.17 


0.1549 


transverse 


\H2 ^ L + 2;H -1 ^ L) (0.4782) 

\H -1 ^ L + 1) (0.4327) 

fii ^ E,H2 ^ L + 8) (0.3867) 


VII 


7.79 


1.2530 


transverse 


|fr - 1 ^ Hi,H2 ^ L + 3)(0.4833) 

\Hi ^ L + 7) (0.3917) 

\Hi L;fi2 ^ L + 8) (0.3791) 


VIII 


8.05 
8.10 


3.5391 
1.1418 


transverse 
transverse 


|fr-2^L)(0.5316) 

\H-1 ^ L + 2) (0.3756) 

\Hi ^ L + 8) (0.3531) 

\H-1^ Hi]H2 ^ L + 3)(0.4788) 

fi2 ^ L + 6) (0.4095) 


IX 


8.87 
8.95 


0.7044 
0.6872 


transverse 
transverse 


|fii ^ L + 11)(0.5061) 

fii ^ E,H2 ^ L + 7) (0.4162) 

\Hi ^ L + 7) (0.4932) 

fi2 ^ E,Hi ^ L + 8) (0.4414) 

fTi ^ L + 4;fr- 1 ^ L + 1)(0.3262) 



GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI2 isomer. 
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Table A. 2: Excitation energies {E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI3 equilateral triangle isomer (c/. Fig. 7), along with the oscillator strength (/12) 
of the transitions. In-plane and transverse polarization corresponds to the absorption due to light polarized in and 
perpendicular to the plane of the triangular isomer respectively. In the wave function, the bracketed numbers are 
the CI coefficients of a given electronic configuration. Symbols H and L, denote HOMO (singly occupied, in this 
case) and LUMO orbitals respectively. HF denotes the Hartree-Fock configuration. 


Peak E (eV) 


/12 


Polarization 


Wave Function 








\HF) (0.8373) 

if- 2 + 5) (0.1329) 


I 3.42 
3.54 


0.0376 
0.1080 


in-plane 
in-plane 


|if-3^L + 5)(0.2908) 
if- 2 ^L + 1) (0.2439) 
if- 2 ^L + 5) (0.3686) 
if-2^ if) (0.3403) 


II 5.61 


0.2565 


in-plane 


|fr-2^L + 5;fr-l^L + 5) (0.4854) 
if^L + l;fr-l^L + 1)(0.4476) 


III 5.87 


0.3413 


transverse 


|fr-3^L + 2)(0.2915) 
if-2^L) (0.2842) 


IV 6.53 
6.53 


6.3289 
5.7925 


in-plane 
in-plane 


|fr^ L + 6) (0.4044) 
fT- 3 ^L + 1) (0.3965) 
fT- 2 ^L + 5) (0.3158) 
fT^ L + 4) (0.3842) 
fr-3^L + 5)(0.2834) 
fr-4^L + l) (0.2256) 


V 6.96 


0.4145 


transverse 


|fr- 2 ^L) (0.3140) 
fT- 3 ^L + 2) (0.2626) 


VI 7.50 
7.57 


0.9430 
0.8630 


in-plane 
in-plane 


|fr-2^L + l;fr^L + 5)(0.3136) 
fr-3^L + 5)(0.2864) 
fr^L + 5;fr-3^L + l) (0.3838) 
fr-3^L + l)(0.2651) 
fr-2^L + 5)(0.2590) 



^ GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI3 equilateral triangle isomer. 
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Table A. 3: Excitation energies (E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI3 isosceles triangle isomer (c/. Fig. 9), along with the oscillator strength (/12) of 
the transitions. In-plane and transverse polarization corresponds to the absorption due to light polarized in and 
perpendicular to the plane of the triangular isomer respectively. In the wave function, the bracketed numbers are 
the CI coefficients of a given electronic configuration. Symbols i^i, H2 and Hs denote SOMOs discussed earlier, H 

and L, denote HOMO and LUMO orbitals respectively. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


GS^ 








\Hl,Hl,Hl) (0.8670) 
if-l^L + 10)(0.1213) 


I 


2.37 


0.0358 


in-plane 


\Hi^ L + l]H^^ L + 2) (0.7066) 
\H -l^ L + l]Hi ^L) (0.4052) 


II 


3.06 


0.0992 


in-plane 


|i^3 ^if2;i^-2^L) (0.4691) 
if - 1 ^ L + l;if3 ^ ifs) (0.4070) 


III 


3.45 


0.0967 


in-plane 


\Hi ^ L + 3) (0.5566) 

\H -l^L + l]Hi ^L) (0.5209) 


IV 


4.11 


0.3208 


in-plane 


\Hi ^ L + 4) (0.6038) 

ifs ^L+l;if-2^L)(0.5272) 


V 


4.83 


0.2242 


in-plane 


\Hi ^L+l;if-2^L+l)(0.5321) 
\Hi ^ L + 5)(0.2611) 


VI 


5.76 
5.85 


5.0792 
0.8553 


in-plane 
in-plane 


|if - 1 ^ L + 1; ifs ^ L) (0.3479) 
if-3^L + l;ifi ^L) (0.2875) 
\H2^ L + l-Hi^ L + 3) (0.2800) 
ifs ^L+l;if-2^L)(0.4081) 
\H -1^ E,H3^ L) (0.2400) 


VII 


5.95 
6.15 


1.7094 
0.7827 


in-plane 
in-plane 


\H -1 ^ L + 2) (0.3296) 

\H -1^L + 1;H3^ L)(0.3138) 

\Hi ^ L + 7) (0.7827) 


VIII 


6.68 


1.7774 


in-plane 


\Hi ^ L + 10) (0.4548) 

\H2 ^ L + l]Hi ^ L + 6) (0.2705) 

\Hi ^ L + 6) (0.2447) 



GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI3 isosceles triangle isomer. 
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Table A. 4: Excitation energies (E) and many-particle wave functions of excited states corresponding to the peaks in 

the linear absorption spectrum of AI3 linear isomer (c/. Fig. 11), along with the oscillator strength (/12) of the 
transitions. Longitudinal and transverse polarization corresponds to the absorption due to light polarized along and 
perpendicular to the axis of the linear isomer respectively. In the wave function, the bracketed numbers are the CI 
coefficients of a given electronic configuration. Symbols i^i, H2 and denote SOMOs discussed earlier, H and L, 
denote HOMO and LUMO orbitals respectively. HF denotes the Hartree-Fock configuration. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


GS" 








\HlH^,H^) (0.8010) 

\H -?> ^ H -l]H -2 ^ L)(0.1913) 


I 


1.24 


0.0317 


longitudinal 


\H2 ^ L + 1) (0.6602) 
i^-l^ifs) (0.3636) 


II 


2.25 


0.0489 


longitudinal 


\H -l^ ifa) (0.6856) 
i7- 2 ^ifi) (0.3230) 


III 


4.01 


0.9019 


longitudinal 


|i7-2^ifi) (0.5249) 
\H -l^ ifs) (0.3471) 


IV 


4.43 
4.47 


2.8593 
0.0960 


longitudinal 
transverse 


\H -l^ ifs) (0.4070) 

\H -l^ L + 4;if2^ iv + 6) (0.2409) 

\H2 ^ L + 2) (0.5402) 

\H -1^ H3]H2^ L + 6)(0.3068) 


V 


4.62 


5.1747 


longitudinal 


li^-l^ifs) (0.4600) 

\H-1^ L + 4]H2 ^ L + 6)(0.2862) 


VI 


5.29 


0.1070 


transverse 


\H2 ^ L + 5) (0.4951) 

\H -1^ H3]H -1^ L + 1) (0.3284) 

if - 1 ^ L + 3)(0.3091) 


VII 


5.83 


0.1412 


longitudinal 


\H -1 ^ L + 2;Hi ^ L) (0.6637) 
if-2^ifi) (0.2225) 
if- l^ifs) (0.2073) 


VIII 


6.31 
6.37 


0.0459 
0.0740 


longitudinal 

transverse 


\Hi ^ L + 6;H3 ^ L) (0.5099) 
\Hi L;if3 ^ L + 6) (0.2706) 
if - 1 ^ L + 3) (0.3989) 
\H -1 ^ H2]H3 ^ L + 6)(0.2266) 


IX 


6.89 


0.1311 


longitudinal 


|ii-5^L + 6) (0.3920) 

if 1 ^ L + 4; ifs ^ i^ + 6) (0.3086) 



^ GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI3 hnear triangle isomer. 
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Table A. 5: Excitation energies (E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI4 rhombus isomer (c/. Fig. 13), along with the oscillator strength (/12) of the 

transitions. In-plane and transverse polarization corresponds to the absorption due to light polarized in and 
perpendicular to the plane of the rhombus isomer respectively. In the wave function, the bracketed numbers are the 
CI coefficients of a given electronic configuration. Symbols i^i,i^2 denote SOMOs discussed earlier, and and L, 

denote HOMO and LUMO orbitals respectively. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


GS* 








\Hl,Hl) (0.8724) 
\H-3^E,H-3^ L)(0.1050) 


I 


1.07 


0.0247 


transverse 


\Hi ^ L + 1) (0.8489) 
if-2^L + 5)(0.1601) 


II 


2.31 


0.3087 


in-plane 


|if-2^ifi) (0.7645) 
\H2 ^ L + 1)(0.3113) 


III 


4.67 


0.5709 


in-plane 


\H -2^ E,H -1^ L + 3)(0.6036) 
i7-l^L + 3) (0.4213) 
\Hi ^ L + 7)(0.3113) 


IV 


4.88 


0.9622 


in-plane 


\H -1^ E,H -1^ L + 3)(0.6036) 
if -3 ^L) (0.4699) 


V 


5.51 


3.8316 


in-plane 


|if-3 ^ L + 4) (0.7378) 
\H -2^ Hi) (0.2161) 


VI 


5.84 


0.4900 


in-plane 


|i7-2^L + 3) (0.3889) 

\H -2^ E,H -3^ L) (0.3758) 

\H-3^ L) (0.3594) 

\H -1^ E,H -1^ L + 3)(0.3591) 


VII 


6.01 


0.5332 


transverse 


\H2 ^ L + 7) (0.7268) 
\H -3^ L + 2)(0.3050) 


VIII 


6.20 


0.7477 


in-plane 


|i7-2^L + 3)(0.5195) 

\H -2^ E,H -3^ L) (0.4189) 


IX 


6.51 


0.2928 


transverse 


|if-3^L + 2)(0.7001) 

\H -2^ Hi]H -1^ L + 2)(0.2232) 

\H -2^ E,H -3^ L + 2)(0.2070) 


X 


6.92 


0.6053 


transverse 


|i7 -3 ^L + 2) (0.5144) 

\H -2^ E,H -3^ L + 2)(0.3549) 

i7-2^L + 5) (0.2676) 


XI 


7.31 


0.4328 


transverse 


|if-2^L + 5) (0.4033) 

\H -3^ E,H -1^ L + 1)(0.3787) 


XII 


7.76 


2.7450 


in-plane 


\Hi ^ L + 8) (0.4387) 

\Hi^ L + 1]H -1 ^ L + 2)(0.3435) 



GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI4 rhombus isomer. 
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Table A. 6: Excitation energies (E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI4 square isomer (c/. Fig. 16), along with the oscillator strength (/12) of the 
transitions. In-plane and transverse polarization corresponds to the absorption due to light polarized in and 
perpendicular to the plane of the rhombus isomer respectively. In the wave function, the bracketed numbers are the 
CI coefficients of a given electronic configuration. Symbols i^i,i^2 denote SOMOs discussed earlier, and iJ, and L, 

denote HOMO and LUMO orbitals respectively. 



r t;cxK 


Hj lev ) 


/12 


Polcirizcition 


Wave Function 










\Hi, Hi) (0.8525) 

\Hi^ L;H -2^ L) (0.0972) 


T 
i 


Z.Uo 


U.UZ I 


in-plcLiie 


\H -l^ Hi]H2^ L + 1) (0.2645) 
\H2 ^ L + 1) (0.2536) 
\Hi L) (0.2443) 


II 


2.68 


0.0301 


in-plane 


\H2 ^ L + 1) (0.4757) 

\ TT -1 > T\(C\ /IQI^Q^ 

\H -1 ^ Hi;H2 ^ L + 1) (0.3608) 


III 


4.19 


0.3420 


in-plane 


|i7-2^L) (0.5889) 

1 TT -1 > T 1 OX/'n /IOQQ^ 

\rL — 1 — > L/ -\- Z) [(j.^Zoo) 
\Hi L) (0.2329) 


IV 


4.92 


0.1131 


in-plane 


\Hi ^ L + 2) (0.5780) 
i/-l^L + 2) (0.4083) 
if-2^L) (0.3198) 


V 


5.33 


U.lzoo 
0.2470 


transverse 
in-plane 


1 zj V T . zj V T \ ^ \ fr\ Q^^OQ^ 

l-O — Z — > L/^ III — > J-J -\- i) [(J.OVlJO ) 

\H -2^ L;Hi^ L + 1) (0.3692) 
\H -2^ Hi;H -2^ L)(0.5193) 
if-l^L + 2)(0.3915) 
i7-2^L + 2) (0.3335) 


VI 


0.00 


1.z44d 


in-plane 


1 TT , r lO\/'^71Q/l^ 

\H -l^ Hi]H -2^ L + 2)(0.2587) 
if- l^L + 2) (0.2579) 


Vii 


D.OO 

6.58 


Q 7QO/1 

0.2634 


in-plane 
transverse 


1 TT , T 1 OX/'n K'7C\(V\ 

\rL — Z ^ L -\- Z)\\j.o(\J\)) 

\H -l^ Hi]H -2^ L + 2)(0.4089) 

\H -I ^ L + 2)(0.3325) 

\Hi ^L + l;if-2^L)(0.4375) 

\Hi ^L + l;if-2^L)(0.4375) 

if- 2 ^L + 3) (0.4183) 


VIII 


6.87 
6.93 


2.9702 
0.2483 


in-plane 
transverse 


|if-2^L + 2)(0.5100) 

if - 1 ^ L + 2) (0.3495) 

ifi ^L + l;if -2^L)(0.3558) 

ifi ^L + l;if -2^L)(0.3558) 

ii-2^L + 3) (0.2929) 


IX 


7.22 


1.4267 


in-plane 


|if - 1 ^ ifi; if - 2 ^ L + 2)(0.4039) 

if -3^ fTi) (0.2900) 

\H -1^ Hi]H -1 ^ L + 2)(0.2848) 



s G5 does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI4 square isomer. 
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Table A. 7: Excitation energies {E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI5 pentagon isomer (c/. Fig. 17), along with the oscillator strength (/12) of the 

transitions. In-plane and transverse polarization corresponds to the absorption due to light polarized in and 
perpendicular to the plane of the pentagon isomer respectively. In the wave function, the bracketed numbers are the 
CI coefficients of a given electronic configuration. Symbols H and L, denote HOMO and LUMO orbitals respectively. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


GS^ 








\iH-2) 
1 1 


^) (0.8679) 


I 


1.03 


0.0195 


in-plane 


\H-l- 
i^ - 1 - 


-> L) (0.8635) 

-> L; if ^ L + 3) (0.0880) 


II 


2.38 


0.0219 


in-plane 


\H-3- 
\H -l- 


->fr- 2) (0.8560) 
-> L + 4)(0.1387) 


III 


3.90 


0.1042 
0.3362 


transverse 
in-plane 


\H ^ L 
\H L; 
\H -4- 
\H -2- 


+ 4) (0.8387) 
H -I ^ L + 2) (0.1944) 
-> L) (0.8140) 
-> L + 9)(0.1841) 


IV 


4.16 


1.3144 


in-plane 


\H -1- 
\H -1- 


-> L + 4) (0.7276) 
-> L + 5) (0.4478) 


V 


4.42 


3.3339 


in-plane 


\H-l- 
\H-l- 
\H -l- 


-> L + 5) (0.7096) 
-> L + 4) (0.4490) 
-> L + 9)(0.1535) 


VI 


4.78 


1.0471 


in-plane 


\H -2- 
\H -2- 


-> L + 9) (0.7992) 

-> L;fr ^ L + 6) (0.2058) 


VII 


5.46 


1.1014 


transverse 


\H ^ L + 13)(0.8156) 

if ^ L; if- 2^ L) (0.1708) 


VIII 


6.37 


0.1270 


in-plane 


|if-3- 


-> L) (0.7632) 


IX 


6.73 


0.7104 


in-plane 


|i/-3- 
\H ^ L 
if - 1 - 


-> L) (0.7370) 
+ 1) (0.3698) 

-> L;fr ^ L + 3)(0.1225) 


X 


7.49 


0.3989 


in-plane 


\H ^ L 
if - 2 - 
\H L; 


+ 3) (0.5087) 
-> L + 16) (0.3508) 
fT-l^L + l) (0.2937) 



GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI5 pentagon isomer. 
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Table A. 8: Excitation energies (E) and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of AI5 pyramid isomer (c/. Fig. 20), along with the oscillator strength (/12) of the 
transitions. In the wave function, the bracketed numbers are the CI coefficients of a given electronic configuration. 
Symbols H and L, denote HOMO and LUMO orbitals respectively. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


GS' 








\{H - 2) ) (0.8591) 

\H -3^ L + 1;H -3^ L + 1)(0.1138) 


I 


1.72 
1.75 


0.0046 
0.0521 


y 

z 


3 ^L + 1) (0.6849) 
i^- 2 ^L + 1) (0.2887) 
\H ^ L + 3) (0.2887) 


II 


2.21 


0.0296 


y 


|if-3^L + l)(0.7170) 
\H - 2 ^ L -\- 2) (0.3402) 
i^- 3 ^L + 2) (0.2290) 


III 


2.55 


0.0477 


z 


\H ^ L-\- 3) (0.5390) 
if-4^if- 2)(0.1296) 


IV 


3.46 
3.48 


0.0399 
0.0769 


y 

z 


\H -3^ E,H -2^ L + 1)(0.6131) 
if- 3 + 2) (0.4975) 
\H - H - 2)(0.7340) 
if- 4 ^L) (0.3735) 


V 


4.04 
4.22 


0.6432 
3.0735 


X 

y 


\H ^ L + 7) (0.5929) 
if ^ L + 4) (0.4432) 
if -3^ L + 2) (0.8272) 
if - 3 ^ L + 1)(0.1580) 


VI 


4.74 


0.3474 


X 


|fr^ L) (0.7617) 
fT^ L + 7) (0.2542) 


VII 


5.08 


0.5494 


z 


\H -2^ E,H ^ L + 2) (0.5540) 
fr-4^L) (0.4833) 


VIII 


5.26 
5.27 


0.3175 
0.1267 


z 

X 


|fr-2^L;fr^L + 5)(0.6251) 
fT- 4 ^L) (0.3902) 
fr-6^fr-2;fr^L + l) (0.6056) 
if L) (0.3242) 


IX 


5.56 


0.1384 


X 


|fr^L + ll)(0.7819) 
fr^L + 13) (0.3051) 


X 


6.00 


1.0052 


X 


|fr^L + 13)(0.8132) 
fr^L + ll)(0.1852) 



^ GS does not correspond to any peak, rather it corresponds to the ground state wavefunction of AI5 pyramid isomer. 
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